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ABSTRACT 

A flight research program is underway at the 
Lewis Research Center to study installation effects 
on turbojet engine exhaust nozzles over a range of 
flight Mach number from 0.60 to 1.30. Two J85-13 
afterburning turbojet engines were mounted in un- 
derwing nacelles. Exhaust nozzle propulsion effi- 
ciency is being measured using load cells, cali- 
brated engines, and flight calibrated nacelle tare 
forces. Reference nozzles of known performance 
were used to determine nacelle tare forces. Analy- 
sis of flight calibration data shows that the sys- 
tem is capable of determining nozzle efficiency to 
a one -sigma randoms rr or of ±1,5 percent. 

INTRODUCTION 

A flight research program is being conducted 
at the Lewis Research Center to study airframe in- 
stallation effects on underwing -nacelle -mounted 
turbojet engine exhaust nozzles at flight Mach 
nimhers from 0.60 to 1.3. This program permits 
large scale testing in the transonic speed range 
where wind tunnels suffer from wall interference 
problems due to blockage effects and shock re- 
flections. The primary objective of the program 
is to evaluate the installed efficiencies of var- 
ious nozzle types for both supersonic cruise and 
supersonic dash aircraft. (l) To do this it is 
necessary to measure accurately engine operating 
parameters so that the internal thrust can be de- 
termined. In addition, a measurement of the noz- 
zle thrust minus drag must be obtained. 

The nozzle internal ideal thrust can be ob- 
tained by using the gas generator method (c.f., 
ref. 2). In this method ground calibrations of the 
engine and afterburner are made. Inflight measure- 
ments of various temperatures and pressure are ob- 
tained and correlated with the ground calibrations 
to obtain the internal nozzle operating conditions. 
An alterna,te approach to determine nozzle internal 
gross thrust is to use a traversing rake at the 
nozzle exit. Results obtained by this method are 
reported in Ref. 3 for a turbojet engine and in 
Ref. 4 for a turbofan engine. In either case ad- 
ditional information is required regarding the flow 
field resiilting from the interaction of the inter- 
nal and external flows in evaluating the overall 
thrust and drag characteristics of a complete ex- 
haust system. For simple exhaust system concepts 
it may be sufficient in determining nozzle drag 
just to measi3re the pressure forces acting upon the 
exterior surfaces of the nozzle boattail. However 
for more complex concepts (such as those utilizing 
auxiliary inlets) many more details of the flow 
influence the nozzle thrust minus drag propulsive 
force. The problem is further complicated when the 
external flow is distorted by the airframe instal- 
lation effects. As a result it may not be practi- 
cal to determine nozzle drag from pressure measure- 
ments since an excessive number of such measure- 
ments would be required. 


The technique to be described in this paper ' 
uses the gas generator method to determine the 
nozzle ideal internal thrust and a load cell to 
measure the nozzle thrust minus drag. An F106 
aircraft (fig. 1) was modified to permit installa- 
tion of two nacelle -mounted J85-GE-13 afterburning 
turbojet engines. The nace3J.es were attached to 
the wings by links which permit forces parallel to 
the nace3JLe axis to be measured by a load cell. A 
calibrated reference nozzle permitted separation 
of drag forces on the foxward part of the nacelle 
from the nozzle forces for the series of flights 
discussed in this paper. The general approach is 
outlined in Fig. 2. Ground calibrations were made 
of the engine, reference nozzle, and load cell. 
These calibrations were used to obtain nacelle 
tare drag forces in flight which finally permitted 
the research nozzle insta3J.ed efficiency to be 
evaluated. 

The purpose of this paper is to describe the 
thrust measurement system, data acquisition sys- 
tem, and to present the results of the calibration 
phase of the program. Special attention will be 
given to the accuracy that can be expected in 
flight using this technique with ground calibra- 
tions of engines and nozzles such as those pre- 
sented in Refs. 5 and 6. An analysis of the ex- 
pected errors was performed and it will be com- 
pared to flight results. 

FLIGHT FACILITY 

Installation 

The underwing nacelles had normal shock in- 
lets and calibrated J85-GE-13 turbojet engines as 
shown in Fig. 3. Also shown on Fig. 3 are the 
front and rear links and the location of the load 
cell. Secondary air to cool the engine and after- 
burner was supplied from the inlet and controlled 
at the periphery of the compressor face by a rota- 
ry valve. 

Thrust Measurement System 

The nacelle support system, consisted of a 
front link, a rear link, and a load ceU. assembly 
located between the two links. The front and rear 
links were each attached to a wing fitting and a 
nacelle fitting with low friction bearings. The 
front and rear links transferred all loads acting 
on the nacelle directly to the wing except loads 
acting in the direction of the nacelle thrust axia 
These loads were transferred to the wing through 
the load cell. 

The load cell was a miniature type containing 
a semiconductor strain-gage bridge and was vented 
so that no ce3J. tare force was encountered due to 
pressure variations. It was attached to a nacelle 
fitting through spherical bushings. 

In order to compensate the load cell for 
forces due to acceleration and nacelle weight com- 
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portents along the nacelle axis, an. accelerometer 
output was combined with the load cell output as 
shown schematically in Fig, 4. The setting of the 
potentiometer determined the portion of the servo 
accelerometer signal which was subtracted from the 
load cell signal to obtain the net thrust minus 
drag output. The potentiometer was adjusted on the 
ground to a position which made this output read 
zero, since the output load cell should be zero at 
static conditions. Low pass filters were incor- 
porated in both the load cell and accelerometer 
circuits to attenuate frequencies higher than the 
natural frequency of the nacelle support system. 

In addition, the thrust measurement was recorded 
48 times during a 11,52 second data scan. . 


would be within ±0.34 percent of full scale. This 
included errors from all sources (ambient tempera- 
ture and press-ure, non-linearity, gain inaccuracy, 
zero offset, drift, noise, etc.) except transducer 
error. To check this accuracy, a reference voltage 
was recorded by the data system during each scan. 
The maximum error of this reading from data re- 
cording d-uring research flights was within ±0.1 
percent of full scale. 

Method of Thrust Measurement 

The parameter of primary interest in the re- 
search nozzle test program is the installed nozzle 
efficiency which is defined as follows: 



To maintain a constant temperature environ- 
ment, the accelerometer mounting block and both 
ends of the load cell were equipped with heaters 
and thermo switches, and the accelerometer and load 
cell were wrapped with insulation. The tempera- 
ture of these units was maintained within 11 K 
(S0° F) or 311 K (100° F) 

In ground calibrations, the accuracy of the 
load cell reading compared with known forces ap- 
plied to the nacelle was within ±0.1 percent of 
full scale. This included hysteresis, nonlinear- 
ity, and data recording system errors, but did not 
include errors due to changes in temperature or 
acceleration compensation. The overall one sigma 
system accuracy is estimated as ±0.25 percent of 
full scale. 


Airborne Data Acquisition System 


The data system was designed to achieve an in- 
herently high accuracy and repeatability. Wherever 
possible, the transducers, instruments, and tech- 
niques used were those which had proven to give 
consistently accurate results with good reliability 
in other programs. Also incorporated in the pro- 
gram was the pre-selection of transducers to obtain 
the best units of a type, and thermostatically- 
controlled electrical heating of the transducers 
was used to minimize' thermal drift during flight. 


The data acquisition system, as shown in 
Fig. 5 consisted of a system to multiplex and re- 
cord quasistatic data in digital form on magnetic 
tape, and a system to reqord dynamic data and var-» 
iations in flight parameters in FM analog form on 
a second magnetic tape. Major components of the 
data system are illustrated in Fig. 6. 


Static pressures were sampled by means of 
scanivalves. The scanivalve unit contained a 
solenoid- actuated rotating pressure passage which 
sequentially connected 48 pressure lines to a 
single transducer. Ten scanivalves located in two 
nacelles provided for a capacity of 480 pressure 
measurements , 


Each of the pressures was recorded only once 
during a data scan. Thermocouples, potentiometers, 
and other transducer outputs were sampled and re- 
corded six times during a scan. The load cell out- 
puts were sampled 48 times during the data scan of 
11.52 seconds. 'The total scan included 1152 words. 


The maximum error specification for the digi- 
tal data system required that 99.73 percent of the 
samples (three sigma, with a normal distribution) 


“ ^)nozzle /T ^ 

T) = (i; 

^ip 

where 

(T - nozzle installed thrust minus drag of the 
nozzle 

F^j^p ideal thrust of the primary jet 

The ideal thrust is obtained from engine cali- 
brations(^) utilizing the gas generator method. 

The airflow is obtained from a typical calibration 
of corrected airflow versus corrected speed, fuel 
flows are metered for both reheat and non-reheat, 
and temperatiu’es and pressures at the discharge of 
the turbine are directly measured. For reheat 
conditions additional calibrations for afterburner 
pressure drop and temperature rise were made. A 
typical plot of afterburner temperature rise vs. 
afterburner fuel-air ratio is shown in Fig. 7. 
Figure 8 shows this inlet, engine, and nozzle com- 
bination in the ground test facility. From these 
conditions and calibrations the ideal primary 
thrust is determined. 

The (T - D) nozzle is obtained from the load 
cell system in the following manner. The load cell 
measures the sum of all forces acting on the na- 
celle. These forces are shown on Fig. 9 and are 
summed as follows: 

+ e. sin e) = (T - 

-(Dram + Dadd. + Dgtrut Dfj^ac Dcowl Dbtmip) (2) 

The summation of terms within the brackets is the 

total of all drag forces (tare force) on the strut 

and nacelle foiward of station 127.68. This is the 

nacelle -nozzle juncture station for all research 

nozzles. The teim mn(a + g sin 9 ) is obtained 

from the output of the accelerometer. When this is 

combined with raw load cell reading, 

left hand side of the equation is referred to as 

the compensated load cell reading. The 

(T - D) nozzle then obtained by knowing the tare 

and the value of the compensated load cell. 

The tare forces have been empirically obtained 
over the range of Mach number, mass -flow ratio, 
and engine power settings by installing reference 
nozzles on each nacelle and flying calibration 
flights. These nozzles were simple, cylindrical 
ejectors with flat external base areas and whose 
internal thrust had been calibrated and whose ex- 
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ternal drag could be easily measured. The internal 
thrust of each reference nozzle was determined by 
correlation with a nozzle thrust coefficient^ Cf . 
Here again it was necessary to use the gas genera- 
tor to establish flow conditions at the engine 
primary exit. Figure 10 shows the functional re- 
lationship of Cf to nozzle pressure ratio for a 
given corrected secondary weight flow of 0.04, for 
various effective primary nozzle areas. The sec- 
ondary air flow was determined from valve open area 
and pressure drop characteristics. The calibra- 
tions were made inflight using ejector nozzle pump- 
ing characteristics determined in a ground cali- 
bration of the nozzle. 

The external drag of the reference nozzles 
was found by pressure measurement on the known 
base area and by calculation of a skin friction 
drag. Skin friction was calculated using flat 
plate theory with the skin friction coefficient de- 
fined as Cf = 0.075/ReV^. 

RESULTS 

The primary objectives of the in-flight cali-. 
bration were to obtain the nacelle tare drags for 
later use in research flights and to assess the 
accTiracy of the thrust measuring system. 

Typical results are presented in Fig. 11. 
Hacelle tare coefficient is plotted versus inlet 
capture mass flow ratio. The data of Fig. 11 were 
obtained at a flight Mach number of 0,90±0.02. 

Data were corrected to Mach 0,90 by using the 
measiared sensitivity of tare coefficient to Mach 
number- in an iterative process. As can be seen on 
Fig. 11 the tare coefficient appears to be a linear 
function of mass flow ratio. This linear relation- 
ship existed for all Mach numbers. A least squares 
curve fit was made of all military power data for 
each nacelle at discrete Mach numbers. Since 
little difference was observed between the left 
and right nacelle tare coefficients they were com- 
bined into one function which was an average of the 
left and right values. 

Cross plots of these data were made at dis- 
crete Mach numbers from 0.60 to 1.30, at a mass- 
flow ratio of 0.82, to obtain tare coefficients as 
a function of Mach number. Fig. 12. The sensitiv- 
ity of tare coefficient to m/mQ was also deter- 
mined and is presented in Fig, 13. Thus the tare 
coefficient can be obtained by determining the 
value at the desired Mach number from Fig. 12 and 
correcting that value to the correct m/m^ by 
using the sensitivity from Fig. 13. 

As can be seen on Fig. 11, the scatter of 
data or deviation from the line becomes greater 
when the engine is afterburning. Only the left 
engine afterburning data are presented but the 
right engine data exhibits the sarpe tendency. The 
afterburning error is evidenced as a bias in tare 
coefficient. That is, all the minimum reheat data 
lies lo-wer than the military data and the maximum 
afterburner data lies above the military. There- 
fore, it is felt that these apparent biases in 
afterburning are the result of errors in the de- 
termination of reference nozzle thrust coefficients 
during calibration and possible inaccuracies in 
the engine afterburner calibrations. An inaccuracy 
in nozzle thrust coefficient in maximum afterburn- 
ing of 1 percent produces a 3.64 percent error in 


tare coefficient. Both of these effects combine 
to produce the apparent shift in drag coefficient 
during afterburner operation. Therefore, these 
biases are not taken into account in determining 
nacelle tare, and the curves of Fig. 11 are based 
on military power data only. 

STATISTICAL RESULTS 

An analysis of flight data was made to deter- 
mine the random error of the tare forces, which 
indicates' the repeatability of the system. The 
random error was evaluated for all military power 
data points "used to determine the tare coefficient 
as a function of m/mQ, This calculation was made 
at various discrete Mach numbers where data were 
available. The results are shown in Fig. 14. The 
standard deviation of the data is defined as fol- 
lows: 



where 

a standard deviation 

Yj, difference between experimental tare co- 
efficient data point and the curve fit of 

data 

n number of data points 

Figure 14 shows the standard deviation as a per- 
cent of tare plotted versus Mach number. Also 
sho-wn as a solid line is the level of error in 
tare necessary to produce a one -percent error in 
nozzle thrust minus drag,' The standard deviation 
of all the data is shown to be within a one per- 
cent error in thrust minus drag. 

The probability distribution of all military 
power data (340 points) over the entire Mach range 
was compared to that of a normal distribution with 
a mean error in drag coefficient of -1.00x10'^ and 
a standard deviation of 20.45x10"'^. This is the 
mean and standard deviation calculated from the 
data. ' If the data were perfectly curve fit by a 
least squares method the mean error would be zero. 
However the curves used were the average of curve 
fits of the left and right nacelles. A plot of 
these distributions is shown in Fig. 15. The 
sample distribution is nearly normal but has a 
greater concentration of data in the region of low 
random error than the normal distribution. The 
sample has 72 percent of the data less than one 
sigma where -as the normal distribution includes 
68 percent of the data less than one sigma. 

ERROR ANALYSIS 

An error analysis was performed to estimate 
the random error in dete mining nacelle drag, and 
nozzle efficiencies. The method used the root -sum 
square of the random errors contributed by each 
measured parameter. 
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where 

E estimated random error for one standard 

deviation 

influence coefficient of any given mea- 
surement 

estimated accuracy of any given measure - 

“7^ ment for one standard deviation 

-^i 

n total nnmber of measurements influencing 

desired re stilt 

Influence coefficients were fotind for each • 
measured parameter using a computer technique. 

Each parameter in the data computation program 
which influenced nacelle drag, or nozzle ideal 
thrust was perturbed a known amotint while all 
other such parameters were held constant, A 
flight data point was used as the base measure- 
ment. The change in nace3J.e tare coefficient, or 
ideal thrust, due to these perturbations was then 
computed. From these changes the influence co- 
efficients were calculated for both military power 
settings and for afterburning. Included in the 
influence coefficient analysis were changes due to 
possible errors in ground calibrations of engine 
airflow, secondary airflow, afterburner pressure 
drop, afterburner temperature rise, and nozzle 
thrust coefficient. 

The estimates of measurement accuracy were 
obtained during calibrations of the various in- 
struments and empirical knowledge of similar mea- 
surements. These estimates were combined with the 
accuracy of the data system to estimate the over- 
all accuracy of the measurement. For determina- 
tion of the ground calibration accuracies the 
basic data were analyzed and the random. error of 
the data were calculated. The one sigma values of 
estimated measurement accuracy are shown in 
table 1. 

Table 2 summarizes the results of the error 
analysis for a military power setting at Mach 
0.90. All parameters affecting nacelle tare force 
and nozzle ideal thiust are listed in the first 
column. They are subdivided into three cate- 
gorieS| measured engine and nozzle parameters, 
calibrated parameters, and measured free stream 
parameters. Influence coefficients which affect 
the tare coefficient are shown in the second col- 
umn, The tare coefficient is the primary end re- 
sult of the present system calibration. The 
third column shows the influence coefficient for 
computed nozzle ideal thrust which is used to 
compute nozzle efficiency. This efficiency is the 
parameter of most interest when the system is used 
to measure performance of a research nozzle. The 
measurement accuracy column was obtained from the 
overall measurement accuracies of table 1 correct- 
ed for the proper range of full scale error. When 
the accuracy was estimated for a calibration 


parameter, such as nozzle thrust coefficient, the 
random scatter in the basic calibration data was 
used. In the final two columns the influence co- 
efficients are combined with the estimated accu- 
racies to obtain the contribution to the final 
tare coefficient and to nozzle ideal thrust. Both 
the tare coefficient and nozzle ideal thrust are 
most affected by compressor face total pressure 
measurement. The random error estimate for the 
tare coefficient was ±3,9 percent and for the noz- 
zle ideal thrust it was ±1.1 percent. 

Tables. 3 and 4 show the results of the error 
analysis for minimum and maximxmi reheat re spec-, 
tively, also at Mach 0.90. For both power settings 
the nozzle ideal thrust is most affected by mea- 
surement errors in compressor face total pressure, 
but the tare coefficient is most affected by pos- 
sible calibration errors in the reference nozzle 
thrust coefficient. This greater error of tare 
coefficient from that at military power setting is 
due in part to the greater thrust level in after- 
burning and also to the calibration accirracy 
assessed to the nozzle thrust coefficient. The 
expected random error in tare coefficient for min- 
imum and maximum reheat was ±5.7 percent and ±6.5 
percent respectively, while the expected random 
error for nozzle ideal thrust was ±1.1 percent and 
±1.2 pel-cent. 

The results of the error analysis of tare co- 
efficient and nozzle ideal thrust were used to 
calculate the expected error in determining the 
efficiency of research nozzles. This calculation 
is summarized in tables 5, 6, and 7. The three 
parameters used to calculate nozzle efficiency 
are tare, load cell force, and nozzle ideal thrust. 
Again influence coefficients were found for these 
three parameters and were combined with accuracy 
values. The accuracy values were taken from the 
tables (E-values). A single accuracy value of 
±3.9 percent for tare coefficient was used for all 
power settings since tare was deterimined by mili- 
tary data only. The expected random error in noz- 
zle efficiencies were calculated to be; ±2.5 per- 
cent for military, ±1.87 percent for minimum re- 
heat, and ±1.81 percent for maximum reheat. 

A coiiparison of random error of flight data 
with, the estimated accuracies of the error analy- 
sis showed that the flight data at Mach 0.90 
(fig. 13) exhibited a random error of ±1.5 percent 
in tare coefficient whereas the estimate was ±3.9 
percent. This difference indicates conservatism 
in the estimates of measurement and calibration 
accuracies. 

Another comparison of maximum afterburning 
data showed a random flight error in tare coeffi- 
cient of ±4.0 percent whereas the estimated error 
was ±6.5 percent. Here too, the estimate showed 
conservatism. 

Since the tare coefficient shows less random 
error in flight than the estimated errors it is 
assumed that the ideal thrust accuracy in flight 
will also be better than estimated. If ;just the 
exhibited flight scatter of ±1.5 percent is used 
in a recalculated nozzle efficiency error estimate 
instead of the original value of ±3.9 percent and 
the same ideal thrust accuracy is used, the error 
estimates are reduced to ±1.5 percent for mili- 
tary, ±1.4 percent for minimum reheat, and ±1.3 
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percent for maxiraum reheat. These estimates are 
probably still conservative because the estimated 
errors in ideal thrust -were not reduced from the 
original estimate. These cal cialat ions are also 
summarized in tables 5, 6, and 7, when the experi- 
mental nacelle parameter is used in the calcula- 
tions rather than -the estimated value, ^ram 


This error analysis was made at only Mach 
0.90 because this Mach number is of prime interest 
in the research program and because more data was 
available at this condition for analysis. The 
error in efficiency should vary somewhat with Mach 
number if the experimental values of tare coeffi- 
cient are used for an error estimate. The error 
should vary according to relationships shown in 
Fig. 13. The error should be greater for Mach 
numbers less than 0.90 and smaller at Mach numbers 
greater than 0.90. 

It should also be recognized that the errors 
quoted are one sigma values so that they are only 
applicable to 68 percent of the data taken. It 
would be necessary to take many data points at 
each condition to be confident that the random 
error of ±1.5 percent in nozzle efficiency would 
be approached. 
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CONCLWIONS 


It was demonstrated that the thrust measuring 
system installed on F106B aircraft is capable of 
determining nacelle tare forces whose repeatabil- 
ity affects nozzle thrust minus drag no greater 
than +1 percent for military power setting. This 
was accomplished by using ground calibrations of P 

engine airflow^, afterburner pressvire drops and 
temperature rises, nozzle thrust coefficient, and p 
a load cell. 

Pq 

An error 'estimate of expected values of noz- 
zle efficiency was made using the experimental 
value of tare coefficient random error of ±1.5 ‘ pg 

percent at a Mach number of 0.90. This estimate 
indicated an expected random error of within ±1,5 
percent in nozzle efficiency for all power q 

settings, for a one sigma deviation, or for 68 
percent of the data. Re . 


SYMBOLS 


S 




primary nozzle exit effective 
area T 


a 

Cp 

Cf 

D 

^add 

^ump 

^cowl 


nace31.e axial acceleration 
ejector thrust coefficient, 

skin friction drag coeffi- 
cient 

nozzle drag 

inlet additive drag 

base drag of ejector 

pressure drag on reward 
facing nacelle step 

pressure drag and skin fric- 
tion drag on inlet cowl 


^5 


Ts 

Vo 

n 

0 


a 


skin friction drag on 
ejector 

skin friction drag on na- 
celle 

ram drag, mV^ 

pressure and skin friction 
drag on strut fairing 

force measured by load cell 

force measured by load cell 
compensated for accelera- 
tion and angle of attack 

acceleration of gravity 

free -stream Mach number 

mass of air captured by in- 
let 

mass of nacelle 

mass of air that could be 
captured by inlet if full 
stream tube were swallow- 
ed. It is based uponfree- 
stream density and ve- 
locity and inlet capture 
area having a diameter of 
37.37 cm (14.715 in.) 

total pressure 

static pressure 

total pressure at primary 
nozzle 

static pressure, internal at 
ejector exit 

free -stream dynamic pressure 

Reynolds number 

reference area 2,2 n? 

(23.9 

ejector internal thrust 

total temperature at turbine 
discharge 

total temperature at primary 
nozzle 

free stream velocity 

nozzle efficiency parameter 
“ ^)nozzle/^ip 

nacelle angle with respect 
to earth horizontal 

standard deviation or band 
within which 68 percent 
of the data falls 
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ratio of secondary total 
temperature at exit of 
secondary passage to pri- 
mary total temperature 

CD ratio of secondary to pri- 

mary -weight flows 
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Type measurement 

Estimated accuracy ~ 1 a 

Engine RPM 

+0.25 percent full scale 

Engine fuel flows 

+0.75 percent full scale 

Temperatures 

1+ 

CJl 

o 

Pressures (scanivalve) 

±0.5 percent full scale 

Altitude 1 

±100 ft 

Mach number 

±0.5 percent 

Free stream total 

±a° F 

temperat-ure 
Load cell 

+0.25 percent full scale 


TABLE 1 
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TABLE 2 MILITAEY POWER, MACH 0.90 
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Figure 3. - Nacelle-engine installation. 
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Figure 4, - Load cell compensation schematic. 








Figure 5. - Data acquisition system schematic. 




























Figure 6. - Data acquisition system components. 






ARERBURNER TEMPERATURE RISE, 


TURBINE DISCHARGE 
TOTAL PRESSURE 

kN/m^abs (PSIA) 

241.2 (35) 

206.9 (30) 

D7.9 (20) 

68.9 (10) 



ARERBURNER FUEL-AIR RATIO 

Figure 7. - J-85 engine after- 
burner calibration. 



Figures. - Propulsion Systems Laboratory's test nacelle. 
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Figure 11. ** Mach number, 0.90. 



Figure 12. - Tare coefficient variation with Mach num- 
ber, at a constant mass flow ratio, 0.82. 
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Figure 13. - Change of take coefficient with mass flow ra- 
tio. 



Figure 14. - Percent random error in tare versus flight 
Mach number, for military power. 
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Figure 15. - Probability distribution of miiitary data. 
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